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SUMMARY 

With the assumption of the linearized pressure-volume relation, 
a solution of the problem of designing an airfoil with any 
theoretically attainable prescribed dimensionless velocity distri- 
bution in a potential flow of a compressible perfect fluid was 
obtained by a method of correspondence between potential flows of 
compressible and incompressible fluids. 

If the prescribed velocity distribution was not theoretically 
realizable, that is, the distribution would result in an open pro- 
file, the method gave an easy way of modifying the velocity distri- 
bution in order to obtain a closed profile. Numerical examples are 
Included. 


INTRODUCTION 

In order to avoid flow separation and excessively high local 
velocities in designing an airfoil or profile, it is advantageous 
to prescribe the velocity distribution as a function of the arc 
length along the airfoil and then compute the airfoil shape. For 
a two-dimensional potential flow of an incompressible nonviscous 
fluid, several methods based on conformal -mapping theory are avail- 
able for obtaining a profile with a prescribed velocity distribution 
in a uniform stream. (See references 1 to 3.) 

A' similar solution for the two-dimensional potential flow of a 
compressible nonviscous fluid was obtained at the NACA Lewis lab- 
oratory. This solution is based on the assumption that the pressure- 
volume relation is given by the tangent to the isentropic curve 
instead of the true curve. The flow pattern of the compressible 
fluid is obtained by transformation from corresponding flow of an 
incompressible fluid in the manner developed by Lin in his f 
extension (reference 4) of the method developed by von Karman and 
Tsien (reference 5) for obtaining the flow of a compressible fluid 
past a given profile. An approach similar to that of reference 4 
is given by Gelbart in reference 6. 
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The method presented herein consists in using the prescribed 
dimensionless velocity distribution about the airfoil in compressible 
flow (l) to select the proper potential flow of an incompressible 
fluid about the unit circle, and (2) to determine the mapping func- 
tion that transforms this incompressible flow into a compressible 
flow. The image of the unit circle under this mapping gives the 
airfoil with the prescribed velocity distribution provided the 
velocity distribution is theoretically attainable. For velocity 
distributions that are not realizable, a method is given for so 
modifying the velocity distribution that a closed profile is 
obtained . 


SYMBOLS 

The following symbols are used herein: 
constants 

complex potential in £ plane 
regular function of J 

\ 

real part of f(£) 
analytic function of £ . 
imaginary part of f([) 
imaginary part ‘ 

regular function of £ 

number determined by included tail angle of airfoil 
P profile 

p pressure 

« 

Q(s) auxiliary function of s 

<1 magnitude of dimensionless velocity in compressible plane 

(ratio of actual velocity to stagnation velocity of sound) 

Re real part 


A,B,C,D, 

a,b 

F(£) 

fU) 

g 

H(£) 

h 

Im 

k 

n 
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8 

w 

w 


arc length on profile 

free -stream velocity In incompressible flow (See fig. 1.) 
complex velocity In incompressible flow 

d 5 


z = x+ly 
a 

r 

7 _ 


& 

i - S+m 
** 

t 

e 

p 

IF 


complex variable (compressible-flow plane) 

angle velocity makes with x-axis (compressible flow) 

N 

circulation (See fig. 1.) , 

angle determined by trailing-edge stagnation point on 
unit circle (See fig. 1.) 

tail angle of airfoil 

complex variable (incompressible-flow plane) 
auxiliary complex variable 
circle angle 
density 

velocity potential 
stream function 


Subscripts! 

c compressible 

1 incompressible 

min minimum 

n nose 

Prime Indicates a related function. 
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THEORY OF METHOD 

/ 


In reference 4, Lin shows that if the pressure -density 
relation for a compressible fluid is expressed as 


P=A-| 

(1) 

where A and B are constants, then given a potential flow of an 
incompressible nonviscous fluid past a profile P^ in the £ 

plane described by the complex potential function F(£ ) and the 
complex velocity w( £) , the potential flow of a nonviscous com- 
pressible fluid about some profile P 0 in the z plane can be 

obtained by choosing a function k( £) , regular and nonzero in the 
exterior of the profile and satisfying the conditions 

•** 

(2) 

on P^ and 

*£ = o 

(3) 

along any path enclosing P^ (The bar denotes the conjugate 
complex quantity.) Then the equations 

z . I + u -/*( t)dt - i/iffy 

(4) 

2< 1 - _ia - 

(5) 

<p e + «r c - »(£) 

(6) 


give the parametric representation of the compressible flov past a 
profile P c in the z plane with [ as a parameter and where 

P c has the same general analytic nature as the original profile 
P i# in these equations cp c , ^ c , q, and a denote for the com- 
pressible flow, the velocity potential, the stream function, the 
magnitude of the dimensionless velocity, and the angle the velocity 
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makes with the x-axis, respectively. Equation (4) gives the 
correspondence between the two planes and equations (5) and (6) 
give the relation between the velocities and the potentials at 
corresponding points. 

It is also shown in reference 4 that if the £ plane is con- 
formally transformed into the *1 plane by the analytical function 
£ * H(£) (taking the profile P^ into a profile P^) and this 
. resulting flow about is mapped into a compressible flow in the 
z plane by 

u *J *( Q 

\_ 

— 2 * e -la = 

1 + /\/i - <1 2 

<P C + i¥ c - F(J) 

where 

P(£) - F(£) 

5(f) - V<f> 

df dj 


M M 

k(£) - k({) 

dj 

then the profile obtained is the same as that given by equations (4) 
to (6). Consequently, the unit circle in the £ plane can be taken 
as P^. When the unit circle is taken as P^, a profile with a 

pointed tall can be obtained as indicated in reference 4 by allowing 
k(£) to vanish at the trailing-edge stagnation point where 

v(£) = 0. 

In order to apply this transformation to the problem of obtain- 
ing the profile with a prescribed dimensionless velocity distribution 
in a compressible potential flow, the complex potential about the 
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circle In the £ plane must he determined and the function k(£) 
obtained. For the actual computation of the profile, only the 
values of these functions for points on the circle are needed. 


Flow in Circle Plane 

The complex potential in the £ plane for the incompressible 
flow about the unit circle £ = e 1® is 

no - q>i + i^i = -w (t + + lo s i + c (7) 

where W is the free-stream velocity and is in the direction of 
the negative £ axis, is the circulation, and C is an 

arbitrary constant. For points on the circle, equation ( 7 ) becomes 

r 

cpi (e) «= -ZW (cos Q + cos 7 ) + — • (0 + it - 7) (8) 

M 

when C is so chosen that ^ = 0 at the trailing-edge stagnation 
point £= -e 1 ?. (See fig. 1 .) The velocity is 

f 

v(e) - -ie“ 10 ^ 2 W sin 0 + ^ ( 9 ) 

Both cp 1 (e) and W( 0 ) are completely determined when I^, W, 

and 7 are known and these quantities are obtained from the pre- 
scribed velocity distribution in the following manner: The magnitude 

of the prescribed dimensionless velocity along the airfoil is given 
as a function of arc length q = q(s) (fig. 2) where, the total arc 
length is taken as 2 a and is measured from the tail along the 
lower surface. If Q(s) is defined as 

Q(s) = -q(s) when 0 s < s n 
Q(s) = q(s) when s n < s < 2 a 
where s n is the nose stagnation point (fig. 2), then 



( 10 ) 
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( 11 ) 


^c^mln 



Q(s) ds 


( 12 ) 


By equation (6) the potentials are equal at corresponding points 
so the 'j 


“ ^c^mlh 


r, =r 

1 c 


(13) 


Then 7 is given by (reference 2) 


cot 7 + y = £ - ”%mln = | _ ^c.mln 

2 Ti 2 r c 


(14) 


and V is given by 


W 


4 n sin 7 4 « sin 7 


(15) 


The flow about the circle is obtained by using these values of 
T^ f w, and 7 in equations (8) and (9) 


Function k( t) 

The function k( £) is computed for points on the unit circle 
by using the prescribed velocity distribution on the airfoil and 
the velocity on the circle to determine the real part of k(£) on 
the circle. Then Poisson’s integral gives the imaginary part of 
k(£) on the circle. The values of k(£) fear any point in the 
exterior of the circle can be obtained from these values on the 
circle by further use of Poisson's integral. For the computation 
of the airfoil, however, only the values on the circle are needed. 

Airfoil with pointed tail. - In order that the airfoil have a 
pointed tail, it is necessary for k(£) to vanish at the trailing- 
edge stagnation point (£ = -e 1 *) of the unit circle. Hence 
k(£) can be written 
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k(t) 


17 \n 


(-f) 


g + ih 


( 16 ) 


where g + ih is a regular function of £ on and outside the 
circle. The exponent n and the included tall angle 5 of the 
airfoil are related by 


_ , 5 

n S 1 • - 

It 


(17) 


as shown in appendix A. For points on the circle 
equation (16) becomes 


£ ■ • i 


k(e) 


=' |l + e 


1(7-6) 


,g(0) + ih(e) 


(18) 


Determination of g(8) and h(9). - From equation (6) the 
potentials cp i (e) and <jf> (s) are equal at corresponding points. 

Thus, by matching these potentials a correspondence is established 
between points along the airfoil arc and the circle angle; that Is, 
s « s(e). By this relation, the magnitude of the dimensionless 
prescribed velocity along the airfoil is obtained as a function of 
the circle angle q = q(e) . Hence, by taking the absolute values 
of equation (5) 


Me) 


1 + 


+ q 2 (e) 


r(e) 

tfe; 


(19) 


for points on the oircle. By use of equations (9) and (18), equa- 
tion (19) reduces to 


2q(e) 


2W sin 0 + 

Clt 


1 + Ji + q2(e) r if 

^ [2 + 2 cob (7 * e)J C 

or, on solving for g(0), 


( 20 ) 


. 8 ( 6 ) 


g(e) a log e (2W sin e + + log € 


A 


1 ± 


2q(e) 


- | l°g 0 + 2cos(7-eJ 

( 21 ) 


The function h(0), the harmonic conjugate of g(e), is 
obtained for points on the circle by using Poisson's integral 
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where D is a constant (reference 7). For convenience, this con- 
stant is taken as zero because a change in D merely rotates the 
coordinates in the z plane, as shown in appendix B. Hence k(£) 
is now completely determined on the circle by using the values of 
g(e) and h(0) (equations (21) and (22), respectively) in equa- 
tion (18). 


Closure conditions. - In order that the airfoil P c be closed, 
k(£) must satisfy the condition 


td)k(t) d£ - » 0 

k(£) 


J 


(23) 


along any path enclosing the unit circle. If these functions are 
expanded about infinity and residues considered, this condition can 
be expressed in a more useful form, which permits the adjustment of 
the values of k(0) as given by equation (18) so that equation (23) 
is satisfied. 


Now 


k(0»(l+Y^ e f(S) 

where f (£) is regular on and outside the circle. By writing 


a i a 2 

f(£) 


P £ f 


then 


k(£) = e 0 


1 + 


ne 


17 + a i / 

1 ■ '■■■ ■ + T>enua xu . 

£ V £J 


+ terms in —r> J > 2 


(24) 


Hence, 


a 

k (£) d £ “ 2n le °( n ® 17 + a q ) 


(25) 
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From the incompressible flow about the circle, 


r (M « fflLfcL » - W 2W1 slny + 

1 d£ l (2 


and 


4 $ ' *' a ° 


W 2 + 


4W^ i sin 7 - mW^e*^ - a^W 


t 


(terms in J >2) 


Therefore 

r\ 


d£ » 2*ie ®(4W 2 i sin 7 - nW 2 e*^ - ajW**) (26) 
Then by using equations (25) and (26), equation (23) becomes 


2sie 


°(ne i7 + a^ - i J^2rtie 0 W 2 (4i sin 7 - ne ly - a^) 


or 


1 Im(ao) e -Eo(ao) ^2M"o> j^i, + + j + 


*p ^4i sin 7 - ne i7 - Re(a 1 ) + i ^(aj^ijy = 0 (27) 

where Re and Im indicate the real and imaginary parts, respec- 
tively. Equating the real part of equation (27) to zero gives 


e 2 -j- n cos 7 + Re(a^) = 0 


(28) 
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4 

But 

Z Ke(a 0 ) _ j£ 
4 


cannot be zero except for zero free-stream velocity; hence 


n cos 7 + Re(a^) = 0 


(29) 


The Imaginary part of equation (27) gives 


2 Re(a Q ) 


n sin 7 + Im(a^) 


y2f 

-1-4 sin 7 + Im(a^) + n sin 7 =0 


or 

Im(a- ) + n sin 7 ^ ■ 8 - ln y (30) 

1 2 Re(ao) . y Z 

T 

Consequently, k( £) vill satisfy equation (23) if Re(a^) and 
IDafa^) satisfy equations (29) and (30), respectively. 

Determination of Re(aQ), Re(a 1 ), and Im^). - In order to 

determine Re(a 0 ), Re(a 1 ), and Im^) when k(0) is given by 
equation (18), the expansion in equation (24) is written as 




/ 
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and 


1 P 2 * 

7j Re k(0) eln 6 d0 = A[n eln y + Imfa^] + b[ii cos 7 + Be^)] 


0 

Then 


jBe(ao) e ^/ A 2 + B 2 


(31) 


A J* Be k ( 0 ) °°8 0 d® + bP 2 * Re k( 0 ) sin 0 d 0 
n cos 7 + Re(a^) * i - Q 0 


.2 t > 2 

A + B 


(32) 


n sin 7 + Im(a. ) ■ — 

A Jt 


r\cn a 2 * 

A Re k( 0 ) sin 6 d£> - B I Re k(e) cob Q d 0 
I ^0 J n 


2 2 
A + B 


(33) 


These equations give Re(a Q ), Be(a x ), and Im^), which are 
required to satisfy equations (29) and (30) for closure. 

Adjustment of k( t) for closure. - If these values of a 0 and 

a-j^ do not satisfy the conditions, k(£) must he modified and this 

modification will change the velocity distribution' slightly In most 
cases. For adjusting k(£), b 1 is defined , 


Im^) 


Re(b^) ■ -jn cos 7 + Re^jJ 


lia^) + n eln 7 W 2 Bin 1 

e 2Re(a 0 ) + y2 
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where Re(a Q ), Re(a 1 ), and Im(a 1 ) are given by equations (31) to 
(33). Then the function k'(£) defined by 

k’(£) - k(£) e (bl/S) (34) 

will satisfy equation (22) and the airfoil obtained using k'(£) 
vlll close. The velocity distribution q' on the resulting profile 
will differ slightly from the prescribed velocity because of this 
adjustment. The ratio of these velocities Is 


al (4|kl 2 -M 2 )*” 8 ^ COB t—ll bl) sln 8 
4 4| k | 2 e 2[Re(bi) cos 9 + Im (bi) sin e] _| v | 2 


Because 


It* (e)l 


kle 


Re(b^) cos 9 + Iin(l)i) sin 0 


it Is possible that k'(£) will not satisfy Inequality (2). If it 
doesn't satisfy the Inequality, a new funtian k”(£) defined by 


k"(£) 


k(J) e 



(35) 


where d Q is chosen so that k”(£) will satisfy Inequality ( 2 ) 
and 


Befd^ = Re(b x ) , r 


Mdi) - Im(b^) - 7 -. * xLsHLJ: 

,2Be(a 0 ) *2 2Be(a + 4 ) £ 

4 4 


will satisfy all requirements and yield a closed airfoil. That such 
a d Q exists is shown by taking d Q to be any positive real number 

such that 


W 2 eln 7 



4 


d 0 >|* e O>i)| + |lm(b^)| 


+ 
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but 


and 


k"(£) 


k(£) 


d + Re(di) cos 0 + 3m(di) sin 0 


|R©(dj_)| +|lm(d^)j<jR©(bi)| +|lm(bi)| 


|Bo(b^ +|lm(b 1 ) 


W 2 sin i W 2 sin y 

2Ro (Sq) y2 IS^dQ y2 


+ T 6 


w 2 


sin y 


£ 

4 


<d. 


Therefore 


k(t)i>ik(t)l>||»(t)| 

which is inequality (2). 

If this method is used, should be chosen as small as 

possible, so that the velocity on the profile will differ only 
slightly from the prescribed values. With k"( £), the Telocity on 
the profile is 


ii 

q = -J 
4 

\iisi 

1 . dQ + Re(d^) cos 6 + Im(dn) sin 9 

|w(e)|e 

k 2 (e)e 2 

dQ + Re(dq) cos 9 + 3m(di) sin 9 

- 

»*<9)| 


In some cases it may be better to modify k'(£) to satisfy 
inequality (2) by changing the higher harmonic terms which, of 
course, do not affect the closure; that is, multiply k' ( £) by 

b 2 bj 

£ 2 ’ £3 ■ ‘ ‘ * 

e .In fact, these terms may be used to reduce 

the changes Induced in the prescribed Telocity distribution when 
modifying k( £) to obtain a closed profile. These terms also may 
be used to alter the airfoil shape when the computed profile is 
too thin or has negative thickness. The selection of the coeffi- 
cients b 2 , b 3 , . . . depends on the particular problem and the 

Telocity distribution desired and therefore no general method oan 
be giTen for determining them. 
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Airfoil Coordinates 


The function k(0) has been obtained to satisfy all require- 
ments; hence, the airfoil coordinates in the z plane are given by 
equation (4) on Integrating around the unit circle. For convenience 
let k(0) be written 


k(0) 


fcj/e) e 


ik 2 (e) 


(36) 


Then from equation (4) the airfoil coordinates are 


z 


y 



(37) 


(38) 


COMPUTATIONAL PROCEDURE 

An outline of the procedure for computing the airfoil follows: 

1. Obtain <P c (s) , T c , and %,mln equations (10), (ll), 

and (12), respectively. 

2. By use of these values in equations (13), (14), and (15), 

obtain 9*, 7 , and W, respectively. Then calculate <£*(0) and 

v(o) from equations (8) and (9). 

3. Plot cp 0 ( B ) and q> i (e) . By comparing equal values of these 

potentials, obtain s as a function of 0, which permits writing 
the prescribed velooity q as a function of 0, q * *(©)• 

4. Comp ute g(0) by equation (21) and h(0) by equation (22). 

Then k(0) is obtained from equation (18). ' 

5. The closure conditions are checked by computing ®®( a g)# 
Re(a^), and Imfa^) from equations (31), (32), and (33), respec- 
tively, Hri ^ checking these values in equations (29) and (30). If 
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these equations are not satisfied, k(0) is adjusted by equation 
(34). This new function k' (©) must he checked in inequality (2). 
If it does not satisfy this condition, k'(0) is further modified 
as indicated in equation (35). 

6. By use of the adjusted values of k(d) that satisfy all 
conditions, the airfoil coordinates are obtained from equations (37) 
and (38) by integration. 


ILLUSTRATIVE EXAMPLES 

Circular profile. - As a check on the method, the theoretical 
velocity distribution on a nearly circular profile, as computed in 
reference 8 from the simplified pressure-density relation, was taken 
as the prescribed velocity distribution (fig. 2) and the profile was 
computed. The profile agrees with the one used in reference 8, as 
shewn in figure 3. 

Airfoil profile. - In this example, the prescribed velocity q 
was obtained from the incompressible velocity distribution about a 
symmetrical Joukowski profile, computed by Lipnan Bers at Syracuse 
University in the form of the ratio of actual velocity to free-stream 
velocity, by taking the dimensionless free-stream velocity to be 0.538. 
The resulting distribution is shown in figure 4 together with the 
final velocity after adjustment for closure. The computed profile 
(fig. 5) is sli^itly thicker than the JoukowBkl profile toward the nose 
and has some reflex camber. The peaks in the velocity distribution 
about the computed profile are lower than those that would occur in a 
compressible flow about the Joukowski profile with the same free-stream 
velocity and angle of attack because: (l) the circulation was kept the 
•same as for the incompressible flow, which resulted in the reflex 
camber; and (2) the thickening of the profile reduced the curvatures 
in the vicinity of the velocity peakB. 


CONCLUSIONS 

By the assumption of the linear pressure-volume relation, a 
method has been given for uniquely obtaining the shape of an air- 
foil having any prescribed velocity distribution, which is theo- 
retically realizable - that is, yields a closed profile - in a 
potential flow of a compressible perfect fluid. If the prescribed 
velocity is not realizable, the method gives a means of modifying 
it so that a closed profile is obtained. Whether the resulting 
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profile is practical or not will depend on other considerations 
such as thickness. The applicability of the method is limited only 
by the accuracy of this linear approximation to the actual pressure- 
volume relation. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 4, 1949. 
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APPENDIX A 

RELATION BETWEEN n AND 6 
On the unit circle, equation (5) becomes 


2q 


4 


-ia e + 2^ (sin 0 + sin y ) 2W 


1 +A7l + q £ 


1 + e 


i(7-e) 


& + ill 


or 


2q 


e 


-ia _ e -1 ^ + 2 ^ (sin 0 + sin y) 2 W 


/ 1 +q 2 

n 

2 i 


y x ~ M 

[2 + 2 cos( 7-0)1 e 

1 + cos( 7 - 0 ) 


+ g 


/ A *1 \ 


For a point on the upper surface of the airfoil very near the 

tail of the airfoil, the angle of the velocity in the z plane is 
given by equation (Al) and 



sin(7-0T) 

+ n arc tan + 

1 + 003(7-0^) 


h(© 1 ) 


where 0 1 is the circle angle corresponding to p^ and 

n <d^<n + 7 

Similarly, for a point p 2 on the lower surface 

3 n sin( 7 - 0 2 ) 

ou = ( 0 ? + — ) + n arc tan 7 ■ - -7 + h( 0 ? ) 

* ■ * 2 1 + cos( 7 - 0 g) c 


where 

- * + ’ ,< 9 2 < - | 
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Hence 


V“s ■ (e r s 2>- * + “ 


arc tao- 


alnCy-gj^) 

1 + COS (7-0^) 


- arc tan- 


sin ( 7-0g ) 

1 + 008(7-0,,) 


h(0 1 ) - h(0 2 ) 


Then the angle 8 Is the limit of (a. -a . ) as p and p 
approach the tall. But 2 


and 


Hence 


limit 

— ► it + 


arc tan 


7 


slnty^) 

1 + 003(7-0^) 


* 

2 


e. 


limit 

+ 7 


arc tan 


sln(7-0 2 ) 


1 + cos(7-0 2 )j 


« 

2 



or 


8 = * + n(-«) 



NACA TN 1913 


21 


APPENDIX B 


INFLUENCE OF D ON PROFILE 


If D is not taken as zero in equation (22), the airfoil 
coordinates x‘ and y' are 


x' = 


^ 2 ( M2 

4k x -^slny^j 
4ki 


sin(k 2 + D + 0) d0 


r\ 2 

4k, 


PtY 

2W sin 7 + r 


4k, 


sin(k + 0) cos D + cos(k 2 + 8)sinD 


d0 


= x cos D - y sin D 


and 


y' =i 


r ,\2 


4kj. - l 2W sin 7 + ^ 


4k, 


cos (kg + D + 0) d0 



+ 0) cos D - sin(k 2 + 0) sin D 


d0 


= y cos D + x sin D 


Consequently, the airfoil given hy x' and y' is the same 
as the one given hy x and y rotated through the angle D. 

,4 ' l 
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